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We report on a numerical reinvestigation of the Aoki phase in lattice QCD with two flavors of 
Wilson fermions where the parity-flavor symmetry is spontaneously broken. For this purpose an 
explicitly symmetry-breaking source term hipijsT^ip was added to the fermion action. The order 
parameter {tpi^^T^'ip) was computed with the hybrid Monte Carlo algorithm at several values of 
(/3, K, h) on lattices of sizes 4* to 12^ and extrapolated to ^ = 0. The existence of a parity-flavor 
breaking phase can be confirmed at /3 = 4.0 and 4.3, while we do not flnd parity-flavor breaking at 
P = 4.6 and 5.0. 



I. INTRODUCTION 

Spontaneous breaking of chiral symmetry is one of the 
main non-perturbative phenomena of QCD explaining 
many features of the hadronic world, in particular of 
hadrons containing u, d and/or s quarks. QCD allows us 
to interpret the light octet mesons as Goldstone bosons. 
The four-dimensional (Euclidean) lattice discretization 
of QCD provides a unique ah initio non-perturbative ap- 
proach. However, in this approach chiral symmetry has 
to be treated with special care. At present, on the lattice 
this symmetry is best realized by satisfying the Ginsparg- 
Wilson relation [1] for the lattice Dirac operator, e.g., 
employing the so-called overlap operator [2, 3], or using 
the five-dimensional domain wall fermion ansatz [4-6]. 
In both cases the Wilson-Dirac operator M^(mo) (with a 
bare mass parameter mo G (—2,0)) serves as an input for 
the fermionic part of the lattice discretized action. 

For the Wilson-Dirac operator (which breaks chiral in- 
variance explicitly) Aoki [7] has argued that in a certain 
range of the hopping parameter k (or the bare mass toq) 
there is a phase in which parity-flavor symmetry is spon- 
taneously broken, in the sense that a condensate as de- 
fined in Eq. (1) exists and is non- vanishing. In agreement 
with the literature we call it the Aoki phase. When k ap- 
proaches the border lines of this phase all pion masses 
tend to zero because one is approaching a second order 
phase transition. In the whole Aoki phase the charged 
pion states are expected to remain massless (in the case 
of Nf ~ 2 flavors) since they appear to be the Gold- 
stone bosons related to parity-flavor breaking, whereas 
the neutral pion should become massive again. The gen- 
eral phase structure as proposed by Aoki is shown in 
Fig. 1. Some numerical results supporting this picture 
were presented in Refs. [7-12]. 

It has been questioned whether the Aoki phase survives 
the continuum limit (in the sense of extending to (3 — oo) 
or, alternatively, ends somewhere at finite /3, perhaps be- 
fore the scaling regime is reached. Previous investigations 
of this problem did not yield a unique answer [10-15]. 



In this paper we present results of a more thorough im- 
merical analysis of this question. As has been discussed 
recently [16], the answer is of relevance for the local- 
ity behavior and the restoration of chiral invariance in 
quenched and full QCD with Ginsparg- Wilson and do- 
main wall fermions. Accordingly, the region of the Aoki 
phase has to be avoided in such computations in order 
not to spoil physical reliability. 

Our investigation was carried out for full lattice QCD 
with Nf = 2 flavors of unimproved Wilson fermions using 
the standard plaquette gauge action. It includes a careful 
extrapolation of {^pij^T^ip) to vanishing external field. It 
shows that the Aoki phase is unlikely to extend beyond 
(3 = 4.6 (which confirms early conclusions in Ref. [14]). 

The outline of our paper is as follows. In Sec. II we 
discuss the proposed phase structure in greater detail. 
Section III provides details of our numerical simulations. 
In Sec. IV we present our numerical results. Section V 
contains the discussion and our conclusions. 



II. THE PROPOSED PHASE STRUCTURE 

Aoki, in his last status report [15], has discussed the 
lattice results supporting the view that for lattice QCD 
with Nf = 2 Wilson fermions there exists a parity-flavor 
breaking phase which is separated from an unbroken 
phase (or from unbroken phases) by second order phase 
transition lines. The conjectured phase structure in the 
{g^, mo) plane is shown on the left-hand side of Fig. 1. As 
can be seen from this figure, two of these critical lines run 
from strong coupling to the weak coupling limit, while 
further critical lines are confined to the weak coupling 
region. At zero coupling, pairs of these transition lines 
join at points referring to the different fermion doublers. 
Aoki [8] has further claimed that along the critical lines 
the pion triplet is massless. The neutral pion becomes 
massless only on the critical lines, due to the presence of 
a second order phase transition, while the charged pions 
turn massless on the critical lines and remain massless 
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FIG. 1: The phase diagram proposed by Aoki et al. in the {g^,Tno) plane (left-hand side) and in the (/3, k) plane (right-hand 
side). The shaded region labeled B denotes the phase where flavor and parity are spontaneously broken. Both symmetries are 
conserved in regions labeled A. 



inside the Aoki phase signahng that flavor symmetry is 
broken. 

When simulating the theory it is natural to draw the 
phase diagram in the (/?, k) plane . Using the well known 
relations k — l/(2mo -I- 8) and (3 = Q/g^, the proposed 
phase structure is mapped to this plane as shown on 
the right-hand side of Fig. 1. Therein the symmetry 
Too —(jno + 8) is hidden in the reflection k ^ ~k 
which is not made explicit for simplicity. The critical 
line Kc{(3) which runs from /3 = to infinity is noth- 
ing but the chiral limit line of lattice QCD. Thus the 
scenario proposed by Aoki et al. might explain why all 
pions are massless along this line despite the fact that 
Wilson fermions explicitly break chiral symmetry. 

In principle, the Aoki phase could be expected to exist 
for all values of (3. In the strong coupling region the exis- 
tence of such a phase was verified by performing numer- 
ical simulations of QCD with Wilson fermions as sum- 
marized in [15] and reconsidered in [12]. For this pur- 
pose a so-called twisted mass term hipij^T^ip was added 
to the action which explicitly breaks parity-flavor sym- 
metry. Without an external field h coupling to t/jij^T^ip, 
the parameter {ipij^T^Tp) would always be zero on a finite 
lattice. {'ipi^^T^'ilj) has to be measured for varying lattice 
size V and non- vanishing h values. The order parameter 
(■!/'i75r^'0)/i=o is then obtained by taking the double limit 
in the following order 

{i>ij5T^i>)h=o = lini hm {tpij^r^Tp) . (1) 

In the literature one finds numerical results from 
quenched [11] and unquenched [10, 12] simulations at 
finite h which support the existence of a parity-flavor 
breaking phase, at least for f3 < 4.0. However, extrapo- 
lations in order to carry out the double limit (1) had not 
been performed. 

Going to larger values of (3 there are contradictory 
statements about the existence of such a broken phase. 
Bitar [14] has come to the conclusion that there is no 
Aoki phase for /3 > 5.0. However, results from quenched 



simulations [11] suggest that the finger structure antici- 
pated by Aoki exists. 

Aoki's scenario was also challenged in Refs. [17-19]. 
In particular, in Ref. [17] it has been argued that flavor 
and parity are not violated at finite lattice spacing. The 
authors have rather proposed that the Aoki phase has 
to be interpreted as a phase with massless quarks and 
spontaneous chiral symmetry breaking. 

In Ref. [13] the controversy has been concisely eluci- 
dated in the sense that, at finite lattice spacing a, the 
Wilson lattice theory is able to exhibit flavor and par- 
ity breaking under certain circumstances. The authors 
have also demonstrated that the results of [18, 19] con- 
cerning the spectrum of the Hermitean Wilson-Dirac op- 
erator 75W(too) (actually obtained for quenched or par- 
tially quenched lattice QCD) lend support, if correctly 
interpreted, to a non-vanishing condensate as defined in 
Eq. (1). 

In terms of an effective chiral Lagrangian, it has been 
pointed out in Ref. [13] that only two possible scenar- 
ios may exist, depending on the sign of one coupling 
coefficient. In the first case, the Aoki picture [8] is ex- 
actly reproduced, whereas in the second case all pion 
masses remain degenerate and non-vanishing over the 
whole ((7^,TOo) plane such that no Aoki phase exists at 
all. If the first case applies to lattice QCD all the way to 
the continuum limit specific predictions concerning the a 
dependence of the neutral pion mass and of the width of 
the Aoki finger pointing towards (toq = 0,g^ = 0) have 
been made. However, if the sign turns to the second case 
the Aoki phase ceases to exist at strong coupling and 
those predictions do not apply. 

After that, the only remaining question is whether the 
Aoki phase really persists until the continuum limit and, 
if it does so, how it shrinks to the point (toq = 0, g^ — 0). 
Only numerical simulations can clarify whether there is a 
strip of parity-fiavor breaking phase in lattice QCD with 
Wilson fermions extending to infinite (3. 
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FIG. 2: Results for (ipi-fsr^^) as a function of k and h aX (3 = 4.0 on a 6* lattice. 



III. SIMULATION DETAILS 



IV. EXTRAPOLATING TO VANISHING 
EXTERNAL FIELD 



We have simulated lattice QCD with two flavors of 
(unimproved) Wilson fermions with (the $ version of) 
the hybrid Monte Carlo algorithm [20, 21] where an 
even/odd decomposition [22] has been employed. An ex- 
plicitly symmetry breaking source term was added to the 
Wilson fermion matrix Af^K? the two- flavor fermion 
matrix was given by 



M (h) = Mw + hi-fr-, 



(2) 



The simulations were performed on lattices ranging from 
4^* to 12^* at (3 values 4.0, 4.3, 4.6, and 5.0, with k and h 
in the intervals 0.15 < k, < 0.28 and 0.003 <h< 0.04, 
respectively. 

In our study we measured {'ipi'y^T'^ ip) & function of 
K at finite h. The parameter (ipi'y^T^ip), which is propor- 
tional to the imaginary part of the trace of 75Af~^(/i), 
was averaged over 100-1000 gauge field configurations 
(see Table II) separated by trajectories of length 1. The 
trace was measured with a stochastic estimator [23]. 

For illustration, results for {ipij^T^iJj) from a 6^ lattice 
at /? = 4.0 are shown in Fig. 2. The location of the peak 
determines the region where subsequent simulations on 
larger lattices and smaller h were performed. In Fig. 2 
the peak is around k = 0.22. It becomes sharper as h 
decreases. We have increased the lattices until measure- 
ments agreed within errors such that we can treat our 
largest lattices as infinitely large. The extrapolation to 
vanishing h is described in the following section. 



In Fig. 3 an analysis of {ipij^T^ip) data is shown for 
/3 = 4.0 and 4.3. As can be seen from the upper and 
lower left plot the interesting region is around k — 0.22 
and K = 0.21, respectively. At these (/3, k) pairs further 
simulations were performed in order to control finite-size 
effects. Data from these simulations are shown in the 
center plots of Fig. 3. No finite size effects are visible in 
the plots except for data from the 4"* lattice at (3 — 4.0. 
Hence, the measurements of (il'ij^T^'ip) from the largest 
lattice at each h can be considered to lie within errors on 
the infinite volume envelope. 

The question arises how to fit these data properly. Mo- 
tivated from the mean field equation 

h ^ Aoa^ + Ai-{k - Kc) a with a = (tpij^T^ip) (3) 



we use the ansatz 



a{h) = A + Bh'^ 



(4) 



It is instructive to look at so-called Fisher plots [24, 25] 
(see the right-hand side of Fig. 3). From Eq. (3) one ex- 
pects data for k < Kc to lie on straight lines ending at the 
origin or at the abscissa, while within the broken phase 
they should lie on straight lines ending at the ordinate. 
As can be seen from the Fisher plots obtained the data 
do not lie on straight lines and therefore do not behave 
mean-field like. 

Using Eq. (4) with the mean field value C — 1/3 re- 
sults in unstable fitting functions, but taking C as a free 
parameter instead, the ansatz describes the data well. 
In fact, the parameter of interest A is robust against 
the introduction of linear and quadratic correction terms 
(see Table I). Furthermore, the fit parameters B and C 
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FIG. 3: In the left column data for {%j)i'ysT^ip) from a 6* lattice are shown as a function of k at several values of h (the lines are 
spline interpolations to guide the eye). The extrapolation to /i = in the infinite volume limit is shown in the center column 
of this figure. The right column shows the Fisher plots with the corresponding fitting function. The upper row shows results 
for P = 4.0, the lower one for /3 = 4.3. 



agree within errors for both values of /3, even when intro- 
ducing corrections. We conclude that the order parame- 
ter {ipij5T^il')h=o is non-zero at (/?, k) = (4.0,0.22) and 
(4.3,0.21). 

Measurements for {ifji^^T^^) at /3 = 4.6 are shown in 
the upper row of Fig. 4. Looking at the upper left plot 
of the figure one sees that (ipi'j^T^tp) still has a peak at 
finite h. The peak becomes narrower and its position 
is shifted from k = 0.1986 to k = 0.1981 as the lattice 
size is increased from 6^ to lO'*. Taking the results from 
the 10^ lattice at k = 0.1981, a fit using Eq. (4) can be 
performed. However, due to low statistics the point at 



h — 0.005 was discarded and therefore some fit parame- 
ters had to be fixed. Using the fit results from the two 
lower values of f3, the parameter A, B and C were al- 
ternately fixed to reasonable values. The extrapolation 
is consistent with a vanishing order parameter (see Ta- 
ble I). The same result is obtained by inspection of the 
Fisher plot in Fig. 4 where the data seem to lie on a 
line ending on the abscissa. This means that the order 
parameter {"ipijsT^ ip) h=o vanishes a,t /3 — 4.6. 

In addition, the parameters B and C agree within er- 
rors for all three values of (3 as can be seen from Ta- 
ble I. Therefore, we also fitted the data globally using 
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FIG. 4: In the left column data for {ipi'y^T'^ip) from lattices of various sizes are shown as a function oi k, ak, h — 0.005 fixed (the 
lines are drawn to guide the eye). The extrapolation to /i = in the infinite volume limit is shown in the center column. In 
the right column the corresponding Fisher plots are shown. The upper row contains data at /3 = 4.6 from lattices of sizes 6^, 
8^ and 10^. The lower row shows measurements for (3 — 5.0 from lattices of size from 6* to 12*. 



ansatz (4) where B and C are common to all data, while 
the parameters Ap and Dp are different for each (3. In Ta- 
ble III the fit results are shown. In agreement with the re- 
sults presented above, the order parameter {ipij5T^ip)h=o 
is found to be finite at /3 = 4.0 and (3 — 4.3, while it van- 
ishes at (3 = 4.6. Furthermore, their values are robust 
against the introduction of a correction term linear in h, 
while B and C are sensitive. 

At /3 = 5.0 a vanishing order parameter becomes man- 
ifest. As can be seen from the lower row of Fig. 4 there 
is still a peak. However, the extrapolation of {ipij^T^ tp) 



to /i = at K = 0.18 as well as the Fisher plot do not 
support a finite value of ('0«75T'^V)/i=o at (3 = 5.0. 



V. DISCUSSION 

In this study we have investigated how far a parity- 
flavor breaking phase in lattice QCD with two flavors of 
dynamical Wilson fermions at zero temperature extends 
in p. An explicitly symmetry breaking term, the twisted 
mass term hipij^T^ip, was added to the Wilson fermion 
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TABLE I: The parameters of the ansatz a{h) =A + Bh^ + Dh-\- Eh'^ fitted to the data for {tpi-y^^r^tp) at /3 = 4.0, 4.3 and 4.6 
with no, linear or quadratic corrections (labeled as 1, 2 or 3). At each h the result from the largest lattice was used in the fit 
(for details see Table II). The data at ft = 0.003 were discarded because these arc from a 6'' lattice. Also the result at 13 = 4.6 
and h = 0.005 was not taken into account due to low statistics. Fixed parameters are presented by their value without giving 
an error. In each case the first fit (bold numbers) was used in Figs. 3 and 4, respectively. 



13 


K 


h 


0.005 


0.01 


0.02 


0.03 


0.04 


4.0 


0.2200 


8^ 


250 


lo-* 


146 


6^ 


1000 






6* 


1000 


4.3 


0.2100 


8^ 


300 




500 




250 


6* 


500 






4.6 


0.1981 






10^ 


170 


10* 


250 


10* 


200 







TABLE II: Statistics used for the final analysis (extrapolation) at selected k values for (3 = 4.0, 4.3 and 4.6. For the respective 
values of h, given in the first row in the second column we report the number of trajectories produced for each lattice size. A 
similar statistic was used for scanning at neighboring /t values. 



fit 


P 




B 


C 


D,/3 


X^/ndp 




4.0 


0.063(2) 











1 


4.3 


0.032(2) 


1.0(1) 


0.64(2) 





3.4 




4.6 


0.004(2) 













4.0 


0.065(2) 






-0.8(3) 




2 


4.3 


0.034(2) 


1.5(2) 


0.71(2) 


-0.8(3) 


2.4 




4.6 


0.004(2) 






-0.7(3) 





TABLE III: The parameters of the ansatz a(h) = Ap + Bh'^ + Dph fitted to the data in Table II for {ipi'ysT'^ip) at /3 = 4.0, 4.3 
and 4.6 with no (fit 1) and linear (fit 2) corrections. The parameter B and C are common to all data, while for each (3 there 
is a separate value for Ap and D/3, respectively. Fixed parameters are presented by their value without giving an error. 



7 



K 





- o 








'<> 


0.23 






inside 


□ : 




o 






border 


O : 


0.22 




o 




vestigial 
















0.21 


o 




o 




- 




o 




\ ^ 






0.20 




o 






- 








o 






0.19 






o 


o 


- 














0.18 








i 












o 





3.6 3.8 4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 

B 
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matrix. The phase diagram was explored in the rectangle 
4.0 < /? < 5.0 and 0.15 <k< 0.28. 

We have presented hybrid Monte Carlo results for 
the order parameter {ipi^^T^il;). The existence of 
a parity-flavor breaking phase could be confirmed at 
(/3, k) = (4.0, 0.22) and (4.3,0.21), where {i^ijsr^i:}, 
measured at finite h, extrapolates to a finite value at 
h — in the infinite volume limit. No parity- flavor break- 
ing was found aX (3 — 4.6 and (3 — 5.0. This suggests 
a phase structure as shown in Fig. 5. Two squares in 
Fig. 5 mark points where we were able to confirm the 
Aoki phase. Two stars mark points where {'ipijsT^iJj) has 
a peak at finite h, but where our extrapolation to /i = 
is consistent with a vanishing order parameter. Conse- 
quently these stars are labeled vestigial. 

According to these results, the Aoki phase for T — 
seems to end close to (3 — 4.6 and k = 0.1981. Rough 
estimates for the upper Kc"^ and lower bound of the 
Aoki phase are 

/3 = 4.0 : 0.215 ~ k^') < 0.220 0.220 < k^") < 0.225 
/3 = 4.3 : 0.205 < < 0.210 0.210 < 4"^ < 0.215 . 

The pair (/3, k) = (4.0,0.215) seems to be quite close to 
the lower boundary. We conclude this from the behav- 
ior of (ipi^^T^ip) in conjunction with the behavior of the 
pion norm [26] (which also has been measured during our 
simulations). At this k) pair {ipi^^T^ip) extrapolates 



to zero at /i = 0, whereas the pion norm seems to diverge 
'AS h ^ 0. Such behavior is expected close to critical lines 

Referring to the discussion of the anticipated phase 
diagram in Sec. II, the results presented here do not in- 
dicate a parity-fiavor breaking phase at /3 > 4.6 which 
was originally claimed to exist at all /3 (see Fig. 1). This 
is suggested not only by the extrapolation of {^i^^T^ijj) 
to /i = in Sec. IV, which yields {i'ij5T^^)h=o = at 
/3 = 4.6 and /3 = 5.0, but also by the observation that 
the peak of {'ipi'-fsT^ip) decreases in height and width as 
f3 increases. The parity-flavor breaking phase seems to 
be pinched off near /3 = 4.6 as illustrated in Fig. 5. From 
the numerical point of view we agree with Bitar [14], who 
has found no evidence of such a broken phase for /3 > 5.0. 

On the other hand, although (ipi^^T^ip) decreases, a 
non- vanishing value at h — in the infinite volume limit 
is not excluded. A decreasing width could have been 
expected from the phase structure in Fig. 1. The fact that 
the peak becomes narrower implies that a high resolution 
scan in k is required at larger values of (3. In addition, 
lattices much larger than 12^ would be needed (which 
is beyond our presently available computing resources). 
With this in mind it is comprehensible that the results 
presented in Ref. [14] could not indicate a broken phase 
for P > 5.0 just because of the small lattice sizes used 
(6"*, S'' and 10^). While we find no numerical evidence 
for the existence of the Aoki phase for f3 > 4.6 one cannot 
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exclude that the phase might be found with methods to 
be invented similar to reweighting. 

A further interesting observation we made is that the 
data behave differently when approaching the parity- 
flavor breaking phase at fixed /3 from k > Kc compared 
with the approach from k < Kc\ First, the peaks of 
{ipi"f5T^ xjj) as a function of k are asymmetric. Second, an 
autocorrelation analysis of {'tpi^^T^ip) shows that mea- 
surements to the right of the peak (above of the Aoki 
phase) are significantly stronger correlated than at all 
smaller k values. 

In light of the possible scenarios discussed by Sharpe 
and Singleton [13] it might be worthwhile to invest more 
computing power in a study of both the width of the Aoki 
finger and the detailed behavior of the neutral pion mass 
inside and outside the Aoki phase with respect to the lat- 
tice spacing dependence. For the case oi Nf = 2 dynami- 
cal (unimproved) Wilson fermions with standard Wilson 
gauge action, however, the impossibility of matching the 
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TT and p masses in the interval 3.5 < /? < 5.3 is known [27], 
which means that scaling is strongly violated. Thus the 
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/3 < 4.6, unfortunately does not allow this potentially 
interesting comparison with chiral perturbation theory. 

In view of the fact that the region above the Aoki phase 
(k > ni^^) is the region of interest for the insertion of 
the Wilson-Dirac operator into the overlap form [3] of 
the massless formion operator an even more extensive 
investigation of this area might be worthwhile to do. 
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